Shear response of a frictional interface to a normal load modulation 
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We study the shear response of a sliding multicontact interface submitted to a harmonically 
modulated normal load, without loss of contact. We measure, at low velocities (V < 100 fim.s -1 ), the 
average value F of the friction force and the amplitude of its first and second harmonic components. 
The excitation frequency (/ = 120 Hz) is chosen much larger than the natural one, associated to the 
dynamical ageing of the interface. We show that: (i) In agreement with the engineering thumb rule, 
even a modest modulation induces a substantial decrease of F. (ii) The Rice-Ruina state and rate 
model, though appropriate to describe the slow frictional dynamics, must be extended when dealing 
with our "high" frequency regime. Namely, the rheology which controls the shear strength must 
explicitly account not only for the plastic response of the adhesive junctions between load-bearing 
asperities, but also for the elastic contribution of the asperities bodies. This "elasto-plastic" friction 
model leads to predictions in excellent quantitative agreement with all our experimental data. 
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I. INTRODUCTION 



Friction between solids carrying a time-dependent normal load is a subject of interest in different 
fields, from mechanical engineering, where the "friction-lowering" effect of external vibrations is 
well known Q and commonly used in applications, to geophysical studies of the effect of rapid 
stress changes on static and dynamic friction of rocks , aiming at a better understanding of the 
coupling between normal and tangential stress states on slipping faults 

These studies involve multicontact interfaces (MCI's), i.e. interfaces between macroscopic solids 
with rough surfaces. The real area of contact thus consists of a large number of small contacts with 
sizes on the micrometer scale. 

In a situation of constant normal load on the MCI, the phenomenological state- and rate-dependent 
friction (SRF) model, formulated by Rice and Ruina jjj] , successfully describes the details of the low- 
velocity dynamics (typically in the 0.01 - 100 /im.s -1 range) of such systems, such as the bifurcation 
between steady sate and stick-slip oscillations M. The model states that the dynamic friction force 
Ff r depends on the instantaneous sliding velocity x and on a dynamic state variable $ as: 



F /r (x, *) = W 



where fio is the dynamic friction coefficient in steady sliding at the reference velocity Vo, and A and 
B are measured to be positive and of typical order 10~ 2 (with B > A). 

The state variable $ can be interpreted as the "age" of the MCI, i.e. the average duration of the 
transient contact between load bearing asperities. For example, in stationary sliding at velocity V, 
the set of microcontacts is destroyed and replaced by a new one over a characteristic sliding length 
Do, and the state variable thus expresses as $ = Do/V. 
More generally, the model specifies the time evolution law of $ as: 

i , X$ 

$ = x - — 

Do 

In Eq. (^) the two corrections to /io have distinct physical meanings: the first term describes an 
instantaneous velocity-strengthening of the interface, while the second expresses strengthening of the 
interface with its "age", which in stationary sliding, where 3> = Do/V, leads to a velocity- weakening 
effect. 

In the spirit of the Bowden and Tabor analysis one can write, for a MCI, the friction force as 
@: 

F /r = o s (x) E r ($, W) 

where a 3 defines an interfacial shear strength, E r is the real area of contact between the solids, 
and W is the normal load carried by the multicontact interface. The age-strengthening effect is 
associated with the creep growth of the microcontact area under normal load: 
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E r ($, W)=Eo(W) 



Eo — the real area of contact at <3> = Dq/Vq — exhibits a linear dependence on W , as explained 
by Greenwood and Williamson's model tVm of contact between rough surfaces. That is, the friction 
force obeys the Amontons law Ff r oc W. 

The velocity-dependent interfacial strength of the interface is described as: 



a a (x) = o s0 



1 + ??ln (l 



This form for the interface "rheology" , discussed in detail in reference |l0| , results from the thermally 
activated depinning of multistable nanometric units localized in a layer of nanometric thickness 
forming a junction between micrometric asperities. 

Eqs. (Q) and yield Eq. ([[]) with /Uo = o s oT,o/W, m = B/fio, r) = A/fio, and since m, r\ <§C 1, 
non-linear logarithmic terms can be neglected. 

The SRF model, and its physical interpretation presented above, have been validated by friction 
experiments on different classes of contacting materials, namely granite, paper, polymer glasses and 
elastomers, under constant normal load applied to the solids. 

In the case of a time-dependent normal load, one can first note that in the Amontons- Coulomb 
description (Ff r = fiW , with constant //), a change in W would lead to a proportional change in 
Ff r , in particular a harmonic normal load modulation W = Wo (1 + e cos (uit)) would produce a 
harmonic frictional modulation about a non- modified average value jiWo- 

In the SRF framework, the variations of x and 3> are non-linearly coupled, through Eqs. (|l|) and 
(Q), to the load modulation, thus resulting in non-trivial effects on the friction force (such as, for 
instance, an anharmonic response to a harmonic normal load) . However, the model as expressed by 
Eqs. (^) and (Q) may not be sufficient to describe correctly the frictional response for the following 
reasons: 

(i) the interface rheology expressed by Eq. (|E|) may not hold for "fast" changes of W, 

(ii) the load variation may modify the interface age strengthening process, thus leading to changes 
in the evolution of E r with $, or in the evolution law of the state variable $ itself. 

Based on their results on the response to normal stress steps and pulses in granite friction ex- 
periments, Linker and Dieterich [0 suggested to modify the evolution law of $, while retaining the 
functional form (^) of the $ dependence of E r . Argueing that a sudden change in normal stress 
would result in a sudden change in $, they postulate: 

$ — 1 - _ — — $, 

D Bo ' 

where they infer a = 0.2, for granite, from their analysis of the response to sudden normal load 
steps. 

In a recent study, Richardson and Marone Q investigated the influence of normal stress modula- 
tions on the so-called "frictional healing" effect in a granular material layer confined between rough 
granite blocks: starting from steady sliding, shear loading is stopped and the subsequent shear stress 
relaxation is measured in presence of a 1 Hz normal load modulation (a modulation frequency close 
to the characteristic stick-slip oscillations frequency that can be inferred for their system). 

Friction experiments with confined granular media have been successfully described by the SRF 
model in situations of constant normal load (although the physical meaning of the variable <3? 
is not clear yet for these systems). However, the use of the constitutive law proposed by Linker and 
Dieterich to include time-dependence of the normal load did not account properly for the details of 
the results obtained by Richardson and Marone. 

In this paper we present an extensive study of the effect of a harmonic modulation of the normal 
load, W — Wo (1 + e cos {tot)), on the dynamic frictional response of a multicontact interface. Ex- 
periments are conducted on an interface between two blocks of poly (methyl-metacry late) (PMMA), 
at velocities V < 100 jj.rn.s~ 1 , load modulation frequency / = 120 Hz ^> V/(2irDo) and relative 
amplitude e in the range 5.10~ 3 — 0.5 (so that no loss of contact between the surfaces occurs). 

We study quantitatively the average F and the components at frequency / and 2/, F\ and F%, 
of the tangential pulling force: 



F = K(Vt — x) = mx + F 



fr 



for different values of V and e ; these results are presented in section [n|. We find in particular 
that the modulation of W induces a systematic decrease of the average dynamic friction coefficient 



p. — F /Wo- This effect, which increases with higher e, is quite substantial: a typical magnitude of 
this effect is a 20% decrease of p, for e = 0.5. 

To analyze quantitatively our experimental data, we need to evaluate which fraction, e e ///e, of 
the load modulation is effectively borne by the microcontacts. Indeed, the normal load modulation 
is too fast for air to be drained in and out of the micrometer-thick interfacial gap. We have studied 
this "leaking air cushion" effect by conducting similar experiments under primary vacuum. From 
these experiments we infer that e e ///e ~ 0.4, and we use this in the subsequent analysis as a scaling 
factor for the modulation amplitude. 

Section is devoted to the analysis of these results in terms of the SRF model and its possible 
extensions to fast load modulations: 

(i) We first test the unmodified SRF model by setting in Eq. ([j]) W — » W{t) and using the 
evolution law (Q) for the state variable $. Numerical integration of these equations leads to a 
quantitatively good prediction of the average friction force F (e). However, the predicted oscillating 
tangential force components Fi and F2 strongly depart from the observed dependences on e and V. 

(ii) We have then tested the proposition of Linker and Dieterich. Using their evolution law ^ 
and their proposed value of a = 0.2, we find that (i) the decrease of F(e) is much smaller than the 
measured values and (ii) the agreement for Fi and F2 is not better than in the previous a = test. 
An attempt with a value of a small enough (a = 0.02) to describe correctly the e dependence of 
F leads to results close to those obtained with the basic SRF equations; this also holds for F\ and 
Fi- That is, as confirmed by a perturbation calculation in e, our experiments are not discriminating 
with respect to the Linker-Dieterich evolution law for such small values of a. 

So, this modified ageing law, even if valid, does not suffice to account properly for the details of the 
frictional response. 

(iii) We propose to modify expression (JE]) for the following physical reason: we know from static 
measurements Jl^] that a MCI exhibits, at shear forces much smaller than the static threshold, 
an elastic tangential response. One can deduce from this a shear stiffness n aS p with the partic- 
ular feature K aap oc W . Now, in our interpretation of friction, the rate variable appearing in o s 
must be the true rate of irreversible (plastic) strain of the interfacial junction of nanometer thick- 
ness h. When taking into account the asperity elasticity n a s P , strictly speaking, this strain rate 
reads h^ 1 d(x — F/K asp ) /dt. In quasi stationary motion, this reduces to the x/h strain rate, hence 
the usual <j 3 (x) expression. In the present experimental situation, K asp is modulated as W itself, 
and the difference between the total and plastic strain rates becomes relevant. Indeed, we show 
that this extended phenomenological elastoplastic generalization of interfacial dissipation leads to 
a very satisfactory description of the average and oscillating shear responses to fast normal load 
modulations. 



II. EXPERIMENTS AND RESULTS 

A. Experimental setup and methods 

The tribometer is composed of a slider of mass M driven along a track through a loading spring of 
stiffness K, one end of which is pulled at constant velocity V, as schematized in the inset of Fig. |l|. 
The slider and track are made of PMMA |wj with nominally flat surfaces lapped with SiC powder 
to a roughness of order 1 /im, thus forming a multicontact interface. 

A detailed drawing of the setup is given in Fig. |lj We impose the velocity V of the loading point, 
in the range 0.1 — 100/xm.s -1 , by means of a translation stage driven by a stepping motor. The 
tangential load is applied on the slider through a leaf spring of stiffness K = 0.2 N./im -1 , which is 
the more compliant part of the system. The dead weight of the slider is 16 N. The average normal 
load Wo can be set in the range 3 — 16 N with the help of a vertical spring attached to a remote 
point itself translated horizontally at the pulling velocity V through a second translation stage, in 
order to prevent any tangential coupling. 

The normal load modulation is achieved by means of a vibration exciter rigidly attached to the 
slider: a harmonic voltage input of given amplitude and frequency / results in a harmonic vertical 
motion of the moving element of the exciter on which an accelerometer is fixed. An acceleration of 
amplitude 7 of this moving element of mass m induces a normal load modulation on the slider of 
amplitude 7717 at frequency /. We thus obtain a normal load W = Wo (1 + e cos(tjt)) with ui = 2nf 
and e = my/Wo in the range 5.10 -3 — 0.5. 

We use the loading leaf spring as a dynamometer by measuring its deflection AX by means of 
an eddy current displacement gauge. The tangential force applied to the slider is thus F = KAX. 
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We measure the average value of the output voltage of the gauge, and use a lock-in amplifier to 
measure the amplitude of the first and second harmonics of this output signal with respect to the 
harmonic excitation signal. We thus characterize the shear force through its average value F and 
its A.C. components at frequency / and 2/, of respective amplitudes F\ and F2. 

The experiments are conducted according to the following protocole: for a fixed set of parameter 
values Wo and V leading to steady sliding when e = 0, we measure p (0) = F (e = 0) / Wo- The 
normal load modulation is then set at amplitude e, while sliding, and shear force measurements yield 
p = F/Wo, pi = |Fi| /Wo and P2 = I-F2I /Wo- The modulation is then switched off and F (e = 0) is 
systematically remeasured before setting a new value of e, in order to check that no drift occurred 
during the measurement. Moreover, we check that for e / the shear force signal does not exhibit 
low-frequency stick-slip oscillations. 

The experimental results reported below have been obtained with an average load Wo — 7 N 
and modulation frequencies / of 120 or 200 Hz, chosen to be away from any mechanical resonance 
frequency of the setup. 



B. Results 



1. Average dynamic friction 



The effect of the normal load modulation on the average tangential force response is to system- 
atically lower the dynamic friction coefficient. The ratio p = F/Wo decreases as the modulation 
amplitude e is increased. The variation Ap (e) = p (e) ~ p(0), plotted on Fig. |^, becomes larger than 
the experimental noise for e > 0.05, and is then quasi linear with e, though it does not extrapolate 
to at e = 0. 

Fig. |^ displays measurements of p (V) for different values of e. It appears that the only effect of 
an increase of the load modulation amplitude is to shift down the p(V) curve, without changing the 
slope dp/d\n(V). Therefore, within experimental accuracy, Ap (e) is velocity-independent. 



2. A.C. components of the force response 



The oscillating force response to a load modulation at frequency / is found to be weakly anhar- 
monic. We characterize it by the amplitudes of the first and second harmonics pi and p2- The ratio 
pi/px lies typically in the range 0.1 — 0.2. 

The reduced first harmonic /ii = |-Fi|/Wo increases monotonically with e and does not show any 
measurable dependence on the driving velocity, as presented on Fig. ^.a where we plot results at 
V — 1 and lOpim.s^ 1 . fii is of order 10 -3 at e — 0.5, i.e. two orders of magnitude lower than the 
average shift Ap. 

The amplitude of the second harmonic in the shear force response also exhibits a monotonic 
increase with the modulation amplitude. Moreover, p2(e) depends significantly on velocity, the 
measured amplitude of this component being lower for smaller V, as presented on Fig. 0.b. 



3. Role of the interfacial air layer 



All the above results correspond to experiments performed at atmospheric pressure. The PMMA 
surfaces in contact are nominally flat over typically £0 = 7 x 7 cm 2 but their roughness implies that 
air is trapped in an interfacial gap of micrometric thickness ho- Any increase in normal load is borne 
in parallel by the microcontacts and by the interfacial air layer. This excess pressure leads the air to 
leak out of the edges of the sample, the rate of flow being limited by the air viscosity. For instance, 
when trying to lift the slider from the track, a strong succion is experienced. One may therefore 
expect that the air layer plays a non negligible role in the interfacial response to load modulation. 

In order to quantify experimentally this "leaking air cushi on" e ffect, we conducted a set of control 
experiments under vacuum. The setup described in section (HA) was placed in a vacuum chamber 
allowing to work at pressures down to 1 mbar (a pressure at which the mean free path of air molecules 
becomes of order 10 pm, i.e. much larger than the interfacial gap, ensuring that the air effect has 
become negligible). 
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We first measure the average dynamic friction coefficient under constant normal load and find 
/Uo ~ 0.5, a value equal to the friction coefficient at atmospheric pressure. This confirms that when 
the interfacial air layer is simply sheared, the corresponding viscous force is negligible with respect 
to the solid friction one. 

Then, following the protocole described in section II A, we measure p,(e), at V = 10/mi.s -1 and 
/ = 200 Hz. We have not been able to use in this control experiment the frequency / = 120 Hz at 
which all other data have been obtained, due to the presence of a spurious mechanical resonance of 
the vacuum chamber close to 120 Hz. 

A comparison of the average friction coefficient variation Afl measured at P — 1 atm and at 
P — 1 mbar is presented on Fig. ^. Note that for a given modulation amplitude, \Ap,\ is larger in 
vacuum than in air. Moreover, when plotted as a function of e e // = e/2.5, the results obtained at 
P — 1 atm are found to collapse on those at P = 1 mbar (see Fig. 

We present in the Appendix a model calculation of the elastohydrodynamic response of the air 
layer. We show that, in all the range of e used in our experiments, the normal response of the 
interface is linear, hence the ratio e e ///e does not depend on e. Moreover, the estimated order of 
magnitude of this parameter at / = 200 Hz is found to be compatible with the above measured 
value. 



III. DISCUSSION AND MODEL 



In this section we analyse our data within the SRF framework. The three parameters A, B and 
-Do involved in the SRF laws are determined experimentally, at constant load W, using the velocity 
dependence of the friction coefficient fl (ln(V)) and the dynamic characteristics of the response close 
to the bifurcation threshold (this method has been described in detail in reference [[To)). We measure 
for our system A = 0.013 ± 0.005, B = 0.026 ± 0.01 and D = 0.4 ± 0.04 fim. All the numerical 
integrations of SRF laws presented below are performed with this set of parameter values. 



A. Rice and Ruina's model 



Before coming to the question of whether or not the Rice-Ruina (RR) equations themselves should 
be modified in the presence of load modulations, it is reasonable to study first which response is 
predicted by the RR model as such. 

Replacing in Eq. (|l|) W by its instantaneous value, the equation of motion of the center of mass 
of the slider reads: 



Mx = K (Vt - x) - Wo (1 + e cos(wt)) 



where x(t) is the instantaneous position of the center of mass of the slider with respect to the track. 
We assume the evolution law of $ (Eq. |^) to be unmodified: 

$ — 1 — — 

Do 



1. The perturbative regime 

Let us first consider the case where e< 1. We linearize Eqs. (^) and (^|) about the steady sliding 
state at velocity V, e = 0: 

<E> st = D /V x st =Vt- Wo/K 

Setting: 

= $ — $ si = 5Re (eXl exp (iuit)) 



8x = x — x st 

we get to first order in e 



= Ke (e$i exp (tut)) 



( o A\ BV 

(-ML? +K + iWoUJyj Xl + Wqjj-®! = -Wop(V) 



i « J f 1+ ( <w + : ^)* 1 = 
where p (V) = no - (B - A) In (V/Vo). We thus obtain 

X 1= ___^ + _ 



1 FA 



where A reads 

KV ( LU\ 2 



A Dn 



and (see reference |1Q|) 



_ (g-A)Wb 
Do 



V 

A D~ 

are, respectively , the c ritical stiffness and pulsation at the stick-slip bifurcation for the unmodulated 
system, uo — y K/M sa 360s" 1 is the inertial frequency. 

In our experimental conditions, with 1/ = 10/im.s -1 and (£> — A)/ A ~ 1, w c — 25s , so w c <C 
u, Wo- On the other hand K/K c > 1. Then A « —uj 2 AWo/V, indicating in particular that inertia 
is negligible. Finally: 

_ KjXij _ epV K 
1X1 ~ Wo 6 ~ ~A^Wo 

Similarly, a second order expansion in e yields the corrections at frequencies 2ui and 0, namely: 

V Kp? 2 
m ~ Lo2W A* e 



Ap ! 



p jp + 2^) 2 



Note that Eq. (|22|) correctly predicts a decrease of the average friction coefficient. 

It is interesting to compare the relative perturbative corrections on the age and velocity variables. 
One finds 



Sx/V 



— < 1 



showing that in our regime, the modulation of the age variable contributes negligibly to the shear 
response. 

Moreover, due to the smallness of A, the effective perturbation parameter is given by: 



A 



that is, the perturbative regime (e< 1(T 2 ) is in practice out of experimental reach. We thus must 
resort to full integration of the above RR equations. 



2. Average friction coefficient decrease 



Numerical results for Afl(e) are plotted on Fig. ^ at excitation frequency / = 120 Hz and 
velocities V = 1 and 10 fj.rn.s~ 1 . Note that a very weak dependence on V is predicted, as observed 
experimentally. 



In section [I B 



we emphasised the role played by the interfacial air layer in our experiments, and 
pointed out that it should be taken into account through an effective modulation amplitude e e //, 
accounting for the fact that only a part of the excitation is borne by the contacting asperities. 
Therefore, the modulation parameter e introduced in Eq. (^) must be understood as e e / / . 

On the other hand, we have measured, at / = 200 Hz, e e f / = 0.4 e. As explicited in the appendix, 
we expect e e / f /'to exhibit some relatively weak frequency dependence. This effect depends crucially 
on the interfacial normal stiffness which is difficult to measure accurately. So, we have chosen to 
treat e e tj /e as a free fitting parameter with an initial trial value 0.4. 

Fig. M shows the best fit obtained for Ap,(e) at V = 1 fj.rn.s~ 1 . It corresponds to e e ff = 0.48 e. 
From now on, all experimental data will be plotted versus this effective modulation parameter. 



3. A.C. response 

The computed first and second harmonics of the frictional response are plotted on Figs, ^.a and 
^.b. One can first notice that the quasi-linear dependence of fji and \xi on V predicted by the 
perturbation calculation also holds here. Moreover, both computed harmonics saturate at large e. 
None of these features agrees with the experimental behaviour. 

We therefore conclude that, in spite of the excellent agreement between the predicted and observed 
Ap (e, V), the unmodified RR model is insufficient to describe the full response of the interface. 



B. Linker and Dieterich's ageing law 

As mentioned in section [j], Linker and Dieterich || (LD) have proposed an extended version of 
the RR model in which the evolution law of the age variable $ is modified according to Eq. We 
now study the predictions of this extended model. 

A perturbation calculation to first order in e, using the equation of motion and the age law 
(^), leads to a first harmonic amplitude: 

e\p,-a\ V K 

fJi = 



A lu W 



This expression points to the fact that the dimensionless LD parameter a acts as a correction to 
the bare dynamic friction coefficient p. LD propose for granite a =0.2 - 0.3, i.e. a sizeable fraction 
of p (~ 0.6 for that material). This leads one to expect that such a value should induce significant 
effects on the predicted shear response. However, we have estimated | 15 | an order of magnitude of 
q for a sparse population of microcontacts (Greenwood interface Jll[) ageing under normal load. 
We have considered the two limits of (i) linear viscoelastic and (ii) fully developped plastic creep, 
using parameters compatible with the measured value of the RR parameter B. Both limits lead to 
the same estimate for a, namely one order of magnitude smaller than the LD value. 

In view of this discrepancy, we have chosen to perform numerical integrations of Eqs. Q and ^ 
for various values of a in the range 0.02 - 0.2. 

The magnitude of the load modulation effect on Ap(e) depends strongly on a, as shown on Fig. 
^.a. Whatever a, Ap(e) remains quasi-independent of V, but for a — 0.2 it is significantly smaller 
than the experimental one. 

Moreover, the dependences of fii and /12 on e and V, shown on Figs. [|b and ^.c, as for the RR 
model, clearly disagree with the experimental results. 

a = 0.02 is found to provide a satisfactory fit for Ap. (e ). How ever, this a value is small enough 
for age effects to become negligible, as noticed in section (IIIA1). We indeed check (Figs, ^.b and 



|.c) that the corresponding fj± and 112 are very close to those obtained from the unmodified RR 
model. 

We are thus led to conclude that: 

(i) The LD evolution law with their proposed value of a does not agree with the experimental results. 

(ii) The Ap data permit to set an upper limit on a without, however, allowing to check the validity 
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of the fuctional form of the LD model. Experiments at much lower frequencies (comparable to the 
stick-slip frequency u> c ) would be needed to answer this question. 



C. Extension of the RR model 



The above analysis suggests that in our "high frequency" regime, where the response is controlled 
by the velocity modulations, it is the "rheological" factor a 3 (x) which should be modified. As 
mentioned in section |, o s describes the plastic dissipation occurring in a junction of nanometer 
thickness between contacting asperities, and the rate involved in a s is a rate of irreversible strain of 
this junction. 

It has been shown that when a multicontact interface is submitted to a shear much smaller 
than the static threshold, its response is elastic. Since the asperity "bodies" (which deform on a 
micrometric thickness, of the order of their diameter) are much more compliant than the nanometer- 
thick elastically pinned adhesive joint [juj, it is their response which controls the interfacial shear 
stiffness K asp . This obeys an extented Amontons law: K asp — W/X, with A a length of order 1 fim 
for our surfaces. 

Sliding amounts to depinning of the adhesive joint which becomes dissipative, while the bodies of 
the asperities retain their elasticity. Therefore, we can schematically represent the sliding interface 
as an elastic element of stiffness K asp , accounting for the bulk elastic strain of the asperities, coupled 
in series to a (frictional) dissipative element (see Fig. |ic| ). When this latter is sheared at velocity 
x p i, the corresponding force is F — f (x p i). 

F ~ fcasp X e l = f\XplJ 

with x e i and x p i respectively the elastic and irreversible displacements. The instantaneous velocity 
of the center of mass of the slider thus reads: 

x = x el + Xpl = d{F/ d ^ p) +r 1 (F) 

and the tangential force on the slider finally reads: 

d 



We therefore express the external force using the same functional form as in Eq. (|l|), but the 
argument of the rate-dependent term becomes x p i. In stationary sliding under constant normal load, 
F/ Kasp is constant, hence the usual dependence on x. In the presence of a load modulation, both 
F and n a3p = Wo (1 + ecos (tot)) /A are modulated, and the elastic strain term becomes significant. 

We present hereafter the results obtained from numerical integration of the corresponding ex- 
tended RR equations. Taking into account the above mentioned fact that in our experimental 
conditions, inertia can be neglected, Eq. (jlh becomes: 



F/W = ^(Vt -x) = fio + Aln 



V & 

I j ■ — ! J r — r I i — /.'Mil 

V dt 



L d 1 --^(Vt- 

ftasp 



+ Bln 



\ D ) 



The parameters A, B and Do are set to their experimentally determined values. The length A 
has been obtained from a quasi-static loading-unloading test J13] at various normal loads. We find 
A = 0.62 ± 0.15 /xm. In view of the relatively large experimental uncertainty on this parameter, 
we have integrated Eqs. ( pg| ) and (Q) with A as a free parameter. The best fit, performed on the 
most sensitive data, namely the Mi( £ e//) ones, is found to correspond to A = 0.7 /im, within the 
experimental uncertainty braket. 

While Afl (e) is found to be only very weakly affected by the rheological correction, this extension 
of the model yields predictions for /ii and /j.2 markedly different from those of both the unmodified 
RR and LD models. Namely, their quasi-linear V-dependences are replaced by much weaker ones. 
On the other hand, neither fj,\ nor \i2 exhibit any longer saturation within the relevant range 
e eff < 0.3. 

As appears from Figs, (fill), the global agreement is now excellent, confirming the validity of the 
extended RR model. 



D. Concluding remarks 



This study leads us to the following conclusions: 

On the one hand, from an engineering point of view, the most spectacular effect of modulating the 
normal load applied to a frictional system is to lower significantly the dynamic friction coefficient. 
This occurs as soon as the modulation is applied, even though its amplitude is low enough to ensure 
permanent interfacial contact between the sliding bodies. 

On the other hand, an important aim of this study was to elucidate the question, relevant to 
seismology, of whether the RR model should be modified to describe the frictional response to fast 
variations of the normal stress. We have shown that, in order to study this, it is essential to measure 
and analyze not only the zero frequency component of the response to an oscillatory load, but also 
its harmonic content. 

In the range of frequencies, much larger than the stick-slip one, that we have studied, the shear 
response is controlled by the velocity modulation, that is by the rate- dependent term of the RR 
constitutive law. However, the quantitative analysis of /ii and ^2 data shows that the relevant 
displacement rate is not, for fast load modulations, the slider velocity, but the rate of plastic defor- 
mation of the adhesive frictional joint. 

This confirms our picture of sliding friction as 2D plasticity prelocalized within a nanometer- 
thick adhesive joint coupled to the bulk of the slider through elastic asperities. 
This enables us to extend correspondingly the expression of the rate-dependent part of the RR state- 
and rate-dependent model. 

The question of the precise effect of a load modulation on interfacial age remains at this stage 
open. Indeed, we have concluded that, at least for our system, this effect is certainly much smaller 
than proposed by Linker and Dieterich. However, precisely for this reason, the "high frequency" 
response is not a good tool for investigating this question. This should be adressed through similar 
experiments at low frequencies, close to the stick-slip one. 



APPENDIX A: ELASTO-HYDRODYNAMIC RESPONSE OF THE INTERFACIAL AIR LAYER 

The aim of this appendix is to establish the equation for the vertical motion of the slider (i.e. 
along the z-direction normal to the interface) and to estimate the relative contributions of the forces 
that are involved. We will, as a result, justify the use of an effective amplitude of modulation of 
the normal load, to account for the fraction of the modulation which is borne by the air cushion 
trapped within the interfacial gap. The order of magnitude of this fraction, referred to as e e ///e in 
the text, and which is the only fitting parameter of our model, is checked independently in a control 
experiment, performed in vacuo, and described in the text. 

The motion of the slider along the z-axis is assumed to be decoupled from the sliding motion along 
x. It is parametrized by the width h of the "gap", i.e. the separation between the average planes 
passing through the rough surfaces of, respectively, the track and the slider. When no modulation is 
superimposed to the bare normal load Wo, the width is ho, a value fixed by the deformation of the 
load bearing asperities which are randomly distributed over the interface of nominal macroscopic 
area Eo. 

a. Elastic response of the multicontact interface. According to Greenwood and Williamson's 
model for multicontact interfaces, the number of load bearing asperities and the real area of contact 
increase linearly with the load. This induces a non-linear dependence on load of the gap thickness. 
Experimentally, it has beeen found that ho — h ~ X z ln(W/Wo), with A z a length, the order of 
magnitude of which is given by the roughness of the surfaces in contact (the standard deviation of 
the surface heights, here 1.3 /xm). In the small amplitude linear regime (Ah <C A z ), the stiffness 
k z = W/X z is a constant and the emlastic restoring force in the z-direction reads: 

Fel „ Wo (±lM (Ai) 

This expression has its exact counterpart for tangential motion, as discussed in the text. Shear 
elasticity involves a length A which is expected to be about 1.7 A z for a Poisson ration v = 0.44. 
Therefore, the measured value A = 0.7 /im yields A z ~ 0.4 /im. 
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b. Elasto-hydrodynamic response of the interfacial air layer. When the gap is e.g. nar- 
rowed, air is compressed until being drained out of the interfacial zone. The resulting pressure 
force on the slider will be denoted F v . For the sake of evaluating F v , we will simplify the problem 
and consider a thin layer of air, of viscosity r\ and density po at atmospheric pressure Po, trapped 
between two perfectly smooth discs of radius R, parallel and distant of h <C R. The relative velocity 
is supposed to vanish at z = and z — h, an assumption which is legitimate if the roughness of the 
surfaces is much smaller than the gap width. Brown and Scholz have reported measurement of the 
gap width between macroscopic ground glass surfaces. At low average pressure corresponding to 
10 _J of the Young modulus of the glass, as encountered in our experiments, they have found that the 
gap width is typically 5 times larger than the rms roughness of the statistically identical surfaces. 
This figure is clearly too small for the "smoothed" model of the interfacial gap to be expected to 
provide a very accurate value of the hydrodynamic force, though it is certainly sufficient to estimate 
its order of magnitude. 

An upper bound for the average pressure excess resulting from the motion of the disk is AP = 
eWo/So, with eWo the amplitude of the normal load modulation and So = irR 2 . The macroscopic 
loading pressure Wo/T,o remains of order 10 mbar in the reported experiments, while e is smaller 
than unity. As a result, AP remains much smaller than the atmospheric pressure Po. However, 
the compressibility of the air layer may be of paramount importance, as suggested by the following 
argument. For infinite plates, no leak occurs at the edge of the gap and the response of the layer, 
trapped under the mean pressure Po, is elastic with an overall stiffness: 

K-air = PoSo/ZlO 

For P = 10 5 Pa, E = 49 cm 2 and ho = 6.5 /im, one finds K a ir = 7.5 10 7 N/m, namely one order of 
ma gnit ude larger than the interfacial stiffness k z originating from the load bearing asperities (see 



Eq. Al) at Wo = 7 N. For finite radius R, edge flow will reduce the amount of air to be compressed in 
order to accomodate the change of gap volume. It is therefore necessary to compute the expression 
of F p by taking account the radial, viscosity controlled, Poiseuille flow which results from the density 
(hence pressure) gradient compatible with mass conservation. 
The continuity equation for the radial flow reads: 

9 , 19/ , x 

_ {ph) + __ {rvph) = 

where v(r) is the mean velocity at radius r (averaged across the gap along the z-direction). The 
pressure field is given by the equation of state of the air at pressures close to Po = 1 atm, which is 
assumed to be: 

P = Po 

P PO 

The set of equations is closed by assuming that the flow is of the Poiseuille type, namely is parabolic 
along the z-direction and varies slowly along the radial direction according to: 

h 2 dP 
V ~ ~Vh]~dr 

As mentioned, the pressure modulation rem ains mu ch smaller than Po, and the gap modulation is 



smaller than ho, hence linearization of Eqs. (A3-A5) is legitimate. One therefore sets P = Po + 5P 



p = po + 5p, and h = ho + Sh, with SP <C Po, Sp <C po and 5h <C ho- Eliminating 5p yields the 
following equation for the pressure field: 

h 2 ld_ ( d(8P) \ _S>_Sh 
1277 r dr \ dr ) P ~ h 

where the dot indicates the partial derivative with respect to time. 

Assuming that the normal elastic stiffness k z of the asperities remains linear, the gap modulation 
res ultin g from the normal load one is harmonic and we therefore seek for a complex solution to 



Eq.A6 of the form: SP = SP exp(iu;t) with Sh = Shexp(iujt). Taking into account the boundary 
condition SP = at r = R and the symetry requirement v = 0, hence dP/dr = 0, at r = 0, one 
obtains: 

6P = ~ P 4 [i-TTtfV 
ho Jo(jR) 
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With Jo the Bessel function of zeroth order and 7 a complex constant given by: 



s/2 V p oh 2 

Integration of the pressure field over the interface yields the complex amplitude of the pressure force: 

nShJ 2 (jR) 



SF P = P -kR' 



h J {"/R) 



with J2 the Bessel function of second order. 

The asymptotic limits deserve comments. For -yR — ► 00, J2/J0 — 1, SF P is real and one recovers 



the purely elastic response (Eq.A2) predicted for large R. It also corresponds to the high frequency 
limit for which the air has no time to leak. For jR — > 0, J2/J0 — > {■yR) 2 /8 = —iuj(3riR 2 )/(2Poho), 
hence SF P is purely imaginary and reduces to a linear viscous damping force which could have been 
derivated by assuming a non compressive Poiseuille flow. For intermediate values of 77?, 8F P has 
both a reactive component, which increases the interfacial stiffness, and a dissipative one. 

c. Prediction for the effective amplitude of load modulation. The slider of mass M oscillates 
in the normal 2-direction under the combined action of the load modulation, the restoring elastic 
force resulting from deformation of the load bearing asperities and compression of the air cushion, 
and the damping force resulting from the air flow. The complex amplitude of modulation of the gap 
width Sh is therefore given by: 

{-Mlu 2 + k z )S\ - SF P = eW 



with SF P given by Eq.AE 



The fraction of the load which is effectively borne by the microcontacts is e e ///e = \K z Sh/(eWo)\. 
It reads: 



nR 2 P X z J 2 (7-R) 



Wo h M7R) 



with u>o = \J k z /M. 

The assumption that e e ff / e does not depend on the amplitude of the modulation relies upon the 
fact that both the elasticity and the viscosity remain linear, namely, as previously discussed, that 
Ah <C A z and ho- This reduces to e e ff < 1, a criterion which is always fulfilled in our experiments. 

For W = 7 N, M = 1.6 kg and A z ~ 0.4^m, uo/(2tt) = 530 Hz. Hence, at 120 Hz, the inertia is 
5.2 10 -2 of the elastic restoring force due to the asperities solely. It is clear from the above analysis 
that a key parameter for evaluating the viscoelastic response of the air is the gap width ho- Taking, as 
discussed previously, the conservative value of 5 times the roughness, namely ho = 6.5 /im, 77 = 10~ 5 
Pa.s, R — 3.9 cm, one computes 7.R = 4(1 — i) and |e e ///e| — 0.24, a value of the same order of 
magnitude than the one, namely 0.48, which is found to provide the best agreement between the 
experimental data and the model prediction for Afi. The role of the interfacial air cushion is further 
confirmed by the control experiment performed in vacuo. At a remaining pressure of 1 mbar, the 
elastic stiffness of the air layer falls two orders of magnitude below the multicontact one. Moreover, 
since the mean free path (at 300 K) of the gas molecules is now of order several 10 fim, i.e. larger 
than the gap width, the viscosity of the layer should vanish. Consequently, the effective amplitude 

vacu 
eff 

An _„2 

the data for A/x(e) performed in the air and in vacuo to collapse on a single curve, as explained 
in the text, one makes use of a scaling ratio which reads explicitely: e^Yf / e eff Uurn — 0.20. The 
experimental value is 0.4. 

The fact that, in both cases, the estimated value is smaller than the observed one by the same 
amount may be possibly attributed to some long wavelength modulations of the gap width ho 
which is likely to remain after the lapping process. Microcontacts may be distributed on patches of 
macroscopic area smaller than Eo , separated by regions of much wider gap in which the air would play 
a negligible role. Typically, a patch radius of 2.5 cm, while keeping the other parameters unchanged, 
would account for the observed value — 0.48 at 120 Hz in the air. This would correspond to 

an effective area of 0.4 Eo, a value still large enough for the microcontacts — the number of which 
does not depend on Eo, according to Greenwood — to remain elastically independent. In addition, 
we have assumed a single degree of freedom for the slider, which is certainly a strong requirement 



is essentially ruled by the slider inertia according to: e^™" ra /e ~ |1 - lu 2 /oj' 2 \ 1 = 1-2 at 200 Hz. 
At atmospheric pressure, keeping the nominal value Eo = 49 cm 2 , £e}// e ~ 0.23. When bringing 
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since the slider is left free to find its own seat on the track. It is clear that a small amount of 
"rocking" would promote the air flow and reduce the cushion effect. 

Finaly, normal load modulation induces a tangential oscillating motion of the slider of ampli- 
tude Ax, hence an air shear flow within the gap. The associated A.C. viscous force on the slider 
ijuiAxT, /ho, which has been neglected in our models, must be compared to the leading term in the 
rate dependent friction force for oscillations about the sliding velocity V, namely AWqljAx/V. The 
ratio of both terms is 77E0V '/ '(Woho) ^ 1CT 7 for V = 100 /im/s, therefore shear viscosity of the layer 
is totally negligible. 
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FIG. 1. Main elements of the experimental setup: Translation stage (Drv); Loading leaf spring (Lsp); Displacement gauge 
(Gg); Vibration exciter (Vb); Weighting spring (Spr); Accelerometer (Acc). The inset is the schematic representation of the 
spring-slider-track dynamical system with control parameters W (normal load), V (driving velocity) and K (external spring 
stiffness) . 

FIG. 2. Variation of the reduced average friction force Aft — (F (e) — F )/Wo versus modulation amplitude e at / = 120 Hz. 
Open symbols correspond to two sets of results at V = 1 fim.s -1 and full symbols to two sets at V = 10 /im.s -1 

FIG. 3. Reduced average friction force p, vs. V for various values of load modulation amplitude. Open circles: e = 0; Open 
triangles: e = 0.2; Full circles: e = 0.35; Full triangles: e = 0.5. 

FIG. 4. Amplitude of the harmonic components of the reduced force response at / = 120 Hz: (a) first harmonic fii(e); (b) 
second harmonic /12(e)- Two sets of data are ploted for each velocity: V = 1 /im.s -1 (open symbols) and V = 10/um.s _1 (full 
symbols) . 

FIG. 5. Reduced average friction force A/x(e), for V = lO/im.s -1 and / = 200 Hz at pressures P = 1 mbar (open circles) 
and P = 1 atm (open triangles). The same set of data at P — 1 atm is also plotted (full triangles) as a function of the scaled 
amplitude e/2.5. 
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FIG. 6. Predictions of the RR model for the average friction decrease Ap,(e), at / = 120 Hz. Lines: RR model for 

V — 1/im.s -1 (full), and V = lO/im.s -1 (dashed). Symbols: raw experimental results at V — 1 fim.s -1 (circles) and 

V — lO/im.s -1 (triangles). The experimental data have been averaged over three different runs. The error bars correpond to 
standard deviations on these runs. 

FIG. 7. Determination of e e // from Afi. Full circles: raw experimental data at V = 1 /jm.s" 1 and / = 120 Hz. Open circles: 
the same set of data plotted as a function of e e // = e/2.1 which provides the best agreement with the RR model prediction 
(line). 

FIG. 8. Comparison between experimental data (symbols) and the RR model predictions (lines) at / = 120 Hz for (a) 
/ii(e e //) and (b) /i2(e e //), at V — l/im.s -1 (full line and triangles), and V = lOpim.s -1 (dashed line and circles). 

FIG. 9. Comparison between experimental data at / = 120 Hz and the LD model predictions for (a) A/Z(e e //); (b) /ii(e e //); 
(c) fJ-2{e e ff)- Lines: LD model for V = lpim.s -1 and a = 0.2 (dotted); V = 10/xm.s -1 , a — 0.2 (dashed); V = 10/itm.s -1 , 
a = 0.02 (full). The predictions for V = 1 fim.s -1 and a = 0.02 are not plotted here because they would be undistinguishable 
from the dotted lines. Symbols : experimental data at V — l^im.s -1 (triangles), and V = lO/im.s -1 (circles). 



FIG. 10. Schematic rheological representation of the frictional interface (see text). 

FIG. 11. Comparison between experimental data at / = 120 Hz (symbols) and the extended RR model predictions (lines) 
for (a) A/x(e e //); (b) /Lti(e e //); (c) fJ-2(e e ff) at V — lpim.s -1 (full lines and triangles) and V = lOpim.s -1 (dashed lines and 
circles). 
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